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Abstract 

In this paper we analyze the double Caldeira-Leggett model: the path integral approach to two 
interacting dissipative harmonic oscillators. Assuming a general form of the interaction between 
the oscillators, we consider two different situations: i) when each oscillator is coupled to its own 
reservoir, and ii) when both oscillators are coupled to a common reservoir. After deriving and 
solving the master equation for each case, we analyze the decoherence process of particular en- 
tanglements in the positional space of both oscillators. To analyze the decoherence mechanism 
we have derived a general decay function for the off-diagonal peaks of the density matrix, which 
applies both to a common and separate reservoirs. We have also identified the expected interaction 
between the two dissipative oscillators induced by their common reservoir. Such reservoir-induced 
interaction, which gives rise to interesting collective damping effects, such as the emergence of 
relaxation- and decoherence-free subspaces, is shown to be blurred by the high-temperature regime 
considered in this study. However, we find that different interactions between the dissipative os- 
cillators, described by rotating or counter-rotating terms, result in different decay rates for the 
interference terms of the density matrix. 
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I. INTRODUCTION 

At the beginning of the 1980s, the work of Zurek Caldeira and Leggett (CL) j^, and 
Zeh and Joos [sl played a decisive role in the understanding of the still unsolved phenomenon 
of quantum measurement; more specifically, the collapse of the wave function and the asso- 
ciated decoherence of superposition states Taking the reservoir into account explicitly 
as a quantum ingredient, and analyzing its effect on the evolution of an initial pure state 
into a statistical mixture, these papers shed light on the shadowy interface between micro- 
scopic and macroscopic domains. Although the wave function collapse remains an obscure 
process, despite striking contributions also dating from the eighties |5[ , much is known today 
about the mechanisms leading to decoherence. In the last few decades we have analyzed this 
phenomenon exhaustively, enabling the proposition of a plethora of protocols to circumvent 
it, ranging from quantum error correction codes QECC 6|] and engineered reservoirs [7| to 
dynamical decoupling and relaxation- and decoherence- free subspaces (R-DFSs) 0,13]. 

More recently, it was demonstrated that entanglement shows scaling behavior in the vicin- 



ity of the transition point 



This connection between the theories of critical phenomena 



and quantum information, together with the search for R-DFSs — which encompasses dis- 
sipative coupled systems — has triggered the study of fundamental quantum processes in 
the domain of many-body physics. Apart from the crucial role played by entanglements in 
the understanding of quantum phase transitions [12], the study of the complex dynamics 
of coherence and decoherence of superposition states in networks of dissipative quantum 
systems has also produced interesting results for quantum information theory [isl, Q ■ 



In particular, in Ref. 



id ] a correlation function was introduced to provide the analytical 



conditions for the existence of R-DFSs in a system of interacting dissipative resonators. This 
correlation function measures the reliability of a decoherence-free subspace. Apart from the 
correlation function, Refs. Q, Q, 1^ put forward, as a conceptual novelty, the need to 
consider distinct reservoirs for distinct quantum systems. 

In this study, we consider a system of two interacting harmonic oscillators, in two situ- 
ations: i) when each one is coupled to its own reservoir and ii) when both oscillators are 
coupled to a common reservoir. As argued in Ref. 10|, the former case, where each system 
interacts with its own reservoir, is the most usual situation. Considering, for example, a 
network of coupled cavities, even when they have the same quality factor their damping 
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mechanisms are independent, except when they interact strongly [10']. By strong interaction 
it is meant that the couphng strength A between the N cavities composing the network must 
satisfy the relation A^A ~ cj, being the natural frequency of each cavity mode in the degen- 
erate case, where Ui = ■ ■ ■ = uj^ = uj. In such a strong coupling limit, the results emerging 
in the case of a common reservoir are completely similar to those for the case where distinct 
reservoirs are assumed. Conversely, in the weak coupling limit, where A^A <^ u, separate 
reservoirs must be assigned to each cavity composing the network. 

Even in particular cases where we could, in principle, assign a common reservoir to differ- 
ent quantum systems, such as a sample of atoms inside the same cavity — the multimodal 
cavity playing the role of a common reservoir — such a reservoir turns out to act as distinct 
reservoirs when the transition frequencies of the atoms are significantly far from each other. 
In this situation, each atom interacts with the reservoir modes around its own frequency 
transition, and the absence of (or small) intersection between the reservoir modes addressed 
by distinct atoms makes the common reservoir act as several distinct ones Evidently, 
when the atomic transition frequencies are not sufficiently apart from each other, the overlap 
between the reservoir modes addressed by each atom indicates that they start to interact 
through their "partially common reservoir". Only in the limiting case, where all the atoms 
in the sample have the same transition frequency — and a full overlap between the reservoir 
modes is achieved — is a "completely common reservoir" accessed by the whole atomic sam- 
ple 1^. At this limit, the interaction between the atoms through their common reservoir 
is maximized, opening up the possibility of an interesting feature arising from interacting 
dissipative systems: the R-DFSs. Therefore, the assessment of such R-DFSs — which may 
become indispensable for the implementation of quantum information — requires a com- 
pletely understanding of the dissipative mechanism for coupled quantum systems: either 
through a common or distinct reservoirs. Apart from the emergence of R-DFSs, the subject 
of collective damping effects has recently produced interesting results, such as the nonad- 
ditivity of decoherence rates observed in a network of dissipative oscillators 0, Q, 14], as 
well as in superconducting qubits jisl. Returning to the atomic samples, we stress that 



the collective damping effects coming from two-atom systems [IS], can be directly identified 
with the nonadditivity of decoherence rates. 

The problem of two coupled harmonic oscillators has already been discussed in the lit- 
erature from various perspectives. We first mention a proposal, based on the possibility of 
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performing a reversible coupling between high-Q cavities, to achieve reversible decoherence 
of a mesoscopic superposition of field states jl?]. A theoretical approach for such an ex- 
perimental proposal is given in Ref. jl8|, where a common reservoir is assumed for both 



cavities. In Ref. 



19| a system of two coupled cavities has been analyzed in which just 



one of the cavities interacts with a reservoir. A master equation is derived for the case of 
strongly coupled cavities and it is shown that the relaxation term is not simply the standard 
one, obtained by neglecting the interaction between the cavities. Finally, in Ref. 13|], each 
cavity is assumed to interact with its own reservoir and a detailed investigation is carried 
out for both regimes of weakly and strongly coupled oscillators. In order that the oscillators 
interact only through their direct coupling and not indirectly through their couplings with 
a common reservoir, it is advisable to assume two distinct reservoirs. It is worth mention- 
ing that a general treatment of a network of coupled dissipative harmonic oscillators has 
recently been presented for any topology — i.e., irrespective of how the oscillators are 
coupled together, the strength of their couplings, and their natural frequencies. As in Refs. 



IG 



131, 



14| , the authors start with a general, more realistic, scenario where each oscillator is 



coupled to its own reservoir, and proceed later to the particular case where all the network 
oscillators are coupled to a common reservoir. 

In the present paper, a general form for the interaction between two dissipative oscillators 
is considered and both situations, of distinct reservoirs and a common one, are analyzed. 
Moreover, instead of the master equation approach, we follow the path integral approach 
adopted by CL in their linear response model [2] , with which considerable progress has been 
made on the subject of quantum dissipation in several areas of physics. In fact, with their 
linear response model, CL accounted for the influence of dissipation on quantum tunneling 



in macroscopic systems 



2l| . Quantum Brownian motion has also been approached through 



the influence-functional method 



22 



23(1 . Moreover, the linear response model in 2| has 



3een applied to many topics in solid state physics, for example the dynamics of polarons 
2J] and a particle coupled to a Luttinger liquid [25 1. 

Long before the CL model, functional integral calculations were used by Feynman to 

and, more recently, they have been 



analyze the problem of polarons in a polar crystal 
applied to the problem of bipolarons [27]. As a direct extension of the Feynman polaron 
model, in the path-integral approach to the bipolaron each electron is harmonically coupled 
to a fictitious heavy particle which replaces the virtual phonon cloud. Each electron also 
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interacts with the fictitious particle of the other electron, apart from the Coulomb repulsion 
between the two. Therefore, the double CL model presented here bears some resemblance to 
the bipolaron problem, with our oscillators (reservoirs) replacing the electrons (fictitious par- 
ticles). By replacing the interaction between the two dissipative oscillators by the Coulomb 
repulsion between the electrons, we end up with a dissipative bipolaron-type model. In fact, 
as we demonstrate here, the oscillators interacts indirectly through their common reservoir 
or even their separate ones. 

Before closing this Introduction, we must mention the recent result in Ref . 28| , where the 
authors point out the effective coupling that is induced between two Brownian noninteracting 
particles by a common reservoir. Such an effective coupling depends on the choice made for 
the spectral function of the reservoir. In the present study, working with an ohmic reservoir 
at the high-temperature limit, we find that this induced effective coupling occurs in the case 
of a common reservoir while, as expected, it is absent in the case of distinct reservoirs. We 
also note that such an effective coupling induced by a common reservoir is also pointed out 



in Refs. [id, ll3|, Q, l20j 



Summarizing, in the present paper we employ the path integral approach to treat a 
network of two interacting dissipative harmonic oscillators. Assuming a general form of 
the interaction between the oscillators, we consider two different situations: i) when each 
oscillator is coupled to its own reservoir, and ii) when both oscillators are coupled to a 
common reservoir. We derive and solve the master equation for each case and, in the latter, 
we identify the reservoir-induced coupling between the oscillators, which arises even when 
the original interaction between them is schwitched off. We verify that such a reservoir- 
induced coupling encompasses both dissipative and diffusive terms which couple together the 
variables of both oscillators. These terms thus account for the energy loss of the oscillators 
through each other, apart from a joint diffusive process. From the solutions we have found 
for the master equations we also compute a general expression for the decay rate of the 
off-diagonal peaks of the density matrix of initial superposition states, which also applies to 
both cases of a common and distinct reservoirs. Finally, considering different interactions 
between the oscillators, we analyze the decoherence process of particular entanglements in 
their positional spaces. 
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II. THE DOUBLE CALDEIRA-LEGGETT MODEL 



The CL problem consists in applying the influence-functional method of Feynman and 
Vernon to study the dissipation-fluctuation dynamics of a quantum system {S) interacting 
with a reservoir (R) modelled by a collection of one-dimensional harmonic oscillators. The 
CL Hamiltonian is given by Hcl = Hs + Hr + Hs/r, where the Hamiltonian of the system, 
represented by a particle of mass m and coordinates q and p, is given by 

= £ + (1) 

the Hamiltonian for the reservoir, consisting of a collection of harmonic oscillators {oOk} of 
coordinates {g^} and {pk}, and masses {nik} (the subscript indicating the kth oscillator 
reservoir) reads 



2 

k 



and the interaction Hamiltonian, linear by hypothesis and defined by the coupling constants 
Cfc, is given by 

Hs/R = q^Ckqk- (3) 

k 

The double CL model considered here consists of two quantum systems Si and S2, of 
masses mi and 1712, coupled through the general form 

^5-1/52 = Aligig2 + Ai2giP2 + >^2iq2Pl + ^22PlP2- (4) 

Two different situations arise, however, in the system-reservoir interactions: i) one that 
seems more appropriate for most physical systems, where each oscillator is coupled to its 
own reservoir, and ii) another, which in practice is rather unusual, where both oscillators 
are coupled to a common reservoir. The Hamiltonians governing the evolutions of the two 
coupled dissipative systems are given, respectively, by 

Hi = Y^ {Hs, + Hs,/R, + Hr^) + Hs,/s,, (5a) 

Hu = J2 {Hs, + Hs,/r) +Hr + Hs,/s, , (5b) 
t 

where £,£',£" = 1,2 from here on. 



In the context of a network of cavities coupled by superconducting waveguides [l7|, l29|, l30| , 



distinct reservoirs must be assumed, in general, for distinguishable cavities, even if they 



exhibit equal quality factors, as long as there are no correlations whatsoever between the 



ing field reaching 



10(1, a sufficiently 



reservoirs. The same applies to distinguishable trapped ions or a trave 
distinguishable optical elements. However, as demonstrated in Ref. 
strong coupling between the cavities or the trapped ions lead to a correlation between the 
reservoirs since, as expected, each particular system of the network start to interact with all 
the reservoirs. There are a few particular situations where a set of quantum systems may 
interact with a common reservoir, such as an atomic sample or distinct fields inside a perfect 
closed cavity. In the former case, different atomic transitions couple with different reservoir 
modes and again, the correlations between these reservoir modes define either a common or 
distinct reservoirs, as will be demonstrated below. In the latter, it has been showed that the 
proximity of the distinct field modes sets the strength of the correlation function between 
the reservoir modes, which governs the emergence of both R-DFSs 1^ . 



A. Distinct Reservoirs 



Starting with the case of distinct reservoirs, the term of the Hamiltonian accounting for 
the oscillators and their interactions with the reservoirs is given by 



2mp 



I 2 2 , 

o + -TT^iklik + ^ikqeqek 



Sl/S2- 



(6) 



The Lagrangian associated with Hamiltonian Hi, which defines the action S = Jq Ci dr of 
the Feynman- Vernon theory, is given by 



^2 



(7) 



where C, the Lagrangian associated with the interaction between the two systems, reads 
C = Aiigig2 + A2i/iiig2 (qi - "^^2 j +Ai2yU22gi (q2 - "^^i ] +^^12 (^1 - ^21^2) (^2 - Ai2gi) 



in which 



(8) 



(9) 



1 — A22^l'^2 

stands for reduced masses which, remarkably, arise exclusively from the piP2 coupling be- 
tween the systems. 
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B. A Common Reservoir 



When a common reservoir is assumed, the term of the Hamiltonian accounting for both 
oscillators and their interaction with the reservoir becomes 



Ha 



The Lagrangian following from Hamiltonian Ha, is given by 

= E f ^ - ^^(*) - E + E ? (^^ - ^lil) - ^ (11) 

€ V k J k 

where £ is as defined in Eq. (IHl). 

Before proceeding to the calculations of the propagator for the double CL model through 
the influence-functional method of Feynman and Vernon, we first diagonalize the Hamilto- 
nian describing both coupled systems 

Hs,+s, = E + ^^(9^)) + Hs,/s,. (12) 

This diagonalization is indispensable to define the strength of the interaction parameters 
{A«/} which result in positive values for the normal-mode frequencies. Otherwise, we could 
lave started from a positive-definite Hamiltonian with a lower bound for the energy spectrum 



3, 



14 



20 



31 



32|. 



III. DIAGONALIZATION OF THE COUPLED SYSTEMS 

Assuming, from here on, that the coupled systems Si and 5*2 are harmonic oscillators of 
frequencies ui and 002, the diagonalized Hamiltonian is described in terms of unitary masses 
and the normal-mode coordinates Qi and Pi, defined in Eqs. f l95p . as 

Hs^^s, = ^ E (P" + ^"Q") ■ (13) 

e 

The original masses have been absorbed by the normal-mode frequencies Qj, obtained 
in Eq. (l88l) . which account for the effective interactions gi, given in Eq. (IHTj) . apart from 
the natural frequencies ui and 002- With the condition that the normal- mode frequencies 
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assume positive values, it is straightforward to show that the relations 

7. + m > \gi\ ^ , 14a 

2 2\ ^ 2 / 2 2\ 

g2\ - Ifi'il ) + {UJ1UJ2) > 2u;it^2 (Ifi'il + |fl'2| ) , (14b) 

must be satisfied. 

The coordinates Qi and P^, associated with the normal-mode frequencies, follow from the 
previous generalized coordinates qi and p^, as described by the transformation (1911) . With 
these normal-mode frequencies and coordinates, the full Hamiltonians Hi and Ha becomes: 



+2W ^ [Re (c^/Qf/) + Im {du'Pe')] CekQek 
V m£u;£ ^ 

where we have used the inverse of the transformation Eq. (19T|) . given by 



(15b) 



= \l [Re (c^/Q^/) + Im {diePe')] , (16a) 



V2^^m^E [Re {du>Pi>) - Im (c^'Q^O] > (16b) 



with the coefficients 



cw = ^ [-Ai, in,) + A2, (r],0] e^"^* , (17a) 



and the functions A^/ (^7^/) given by Eq. fIDObp 



Therefore, through the diagonalized form of the Hamiltonian Hsj^+S2> we get two inde- 
pendent and generalized CL Hamiltonians in terms of normal coordinates Qi and Pg. The 
Hamiltonians fllSaj) and (]15b|) describe two independent harmonic oscillators, both interact- 



ing with a common reservoir in the latter case and each one interacting with its own reservoir 
in the former. By a generalized CL Hamiltonian we mean that the system-reservoir coupling 



exhibits, due to the general form of the interaction between the two oscillators in Eq. ([5]), a 
momentum-position term apart from the usual position-position one. Thus, differently from 
the CL model, where^a momentum-position coupling is easily handled through a coordi- 
nate transformation [2], in the present double CL model this momentum-position coupling 
hampers the application of the influence-functional method, since the derived Lagrangian in- 
volves tricky reservoir-reservoir interaction terms. Therefore, in spite of the diagonalization 
of the interaction described by Hamiltonian Hsj^+S2^ is preferable to approach the problem, 
for both distinct reservoirs and a common one, through the Lagrangians in Eqs. ([7]) and 
( ITTj) . respectively. As demonstrated below, in the former case we end up with a product of 
two influence functionals, identical to that of the CL model, whereas for a common reservoir 
the influence functional does not factorize into individual functionals. 



IV. FEYNMAN- VERNON THEORY 



To obtain the reduced master equation describing the time evolution of the coupled 
systems 5*1 and 5*2, we proceed from the integral form of the density operator of the whole 
system at a time t, given by p(t) = U(t)p{0)W{t), where the evolution operator U(t) follows 
from Hamiltonians Hi and Ha. Next, we assign the coordinates Xi and yi to system Si and 
the A^-component vectors = {Xn, . . . , X^^) and Yi = (Y^i, . . . , Y^tv) to reservoirs Ri. 
In the case of a common reservoir, we only have to disregard one of the reservoirs Re, in 
the expression derived for the case of distinct reservoirs, to obtain the associated influence 
functional. Using the notation {xe} = xi,X2 and {dx^} = dxi dx2, and the same form for 
other variables, we obtain the matrix element in the coordinate representation 

({a;,},{X,}|p(t)| W,{Ya) 

= |dK}d{i/;}d{x;}d{Y;} (K},{x;}|p(o)|K},{y;}) 

X K {{xe} , {X,} , t; {x'i} , {Xa , 0) K* {{ye} , {Y,} , t; {y',} , {Y^} , 0) , (18) 

where functional integrations are evaluated over paths xe{t'), ye{t'), ^e{t') and Ye{t'), with 
endpoints xe{t) = xe, xe{0) = x'^, ye{t) = ye, ^^(O) = y'(,, Xe{t) = Xe, X^(0) = X^, Ye{t) = Ye, 
and Ye{0) = Y^. The propagator K is given by 

K {{xe} , {Xe} , t; {x[} , {X',} , 0) = Jd {xe} D {X,} exp (^^S [{xe} , {X,}]^ . (19a) 

10 



The action S follows from the Lagrangian C = oi Ca, a.s S = / Cdt'. Tracing out 

Jo 

the reservoir coordinates, we obtain the reduced density operator describing the coupled 
systems 5*1 and 5*2 under the influence of their respective reservoir, given by 

p{{xe},{ye},t) = Jd{4}d{y',} J {{xe} , {ye} ,t; {x',} , {y',} ,0)p{{x',} , {y',} ,0) , (20) 

where we have assumed that the system-reservoir coupling is turned on suddenly, such that 
the total density operator is initially given by p(0) = p{0)prj^+r^{0). The propagator for the 
density operator turns out to be 



Hi^e} , {ye} ^t; {x'g} , {y'^} ,0) 

D{xe}B{ye} [{xe} , {ye}] 



X exp — 



Yl ^"^s, [{xe}] - Ss, [{ye}]) + Ss,/S2 [{xi}] - Ss,/S2 [{Ve}] 



(21) 



Although expression (!20!) ap plies to both cases, i) and ii), the computation of the influence 
functional JF [{xe} , {ye}] |33|, [SJ] representing the effects of the reservoirs on the systems is 
completely different for the two cases. 



A. Distinct Reservoirs 



When considering distinct reservoirs, the influence functional in Eq. (1211) simply factorizes 



as 



J^[{xe},{yi}] = ll^{^e,yi], (22) 

e 

and the component arising from the interaction of system Se with reservoir Re, is written as 

J'[xe,ye] = J dX^dX^dY^p^, (X;,y;,0) J DX.DY, 

X exp I ^ [X,] - Sr^ [Ye] + Ss,,r, [xe, X,] - Ss,,r, [ye, Y,]) | (23) 

Evidently, when disregarding the interaction between the systems, described by 
Hamiltonian Hs^/S2 Eq. (jl]), the propagator J in Eq. fl^ reduces to a product of 
propagators identical to that obtained in the CL model. However, the influence functional 
factorizes independently of any requirement for non-interacting systems, since it only takes 
into account the interaction between the systems and their respective reservoirs. Therefore, 
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the influence functional (l22l) is obtained directly as a product of the functional obtained from 
the CL model. All that remains to be done is to obtain the propagator fl2Tl) by considering 
the actions coming from the interacting systems. 

Thus, from the CL model we obtain directly the form 

J^[xe,yi] = exp dr^ dt' ae^ir - f) [xe{r) - yi{T)][xi{t') - y(>{t')] 

X exp (-^ / dr [ dt' an{r - t') [a;,(r) - y,(r)] [xi{t') + yi{t')] \ , (24) 







h 

with the real and imaginary parts of a function a (t) given by 



^(^) = cos [uekit)] , (25a) 



k 

«n(t) = - V 7^^^ Sin [.;,,(t)] . (25b) 

k 

Assuming that the reservoir modes are sufficiently closely spaced to allow a continuum 
summation, we define the spectral functions [2I 

The introduction of a frequency cutoff fi^ considerably higher than the characteristic fre- 
quencies of the problem, together with the assumption of an Ohmic reservoir, where the 
distributions XiioS) are defined by the damping constants r^^, such that 

, . J ^ < fif , 
XH^) = < ^ (27) 

I 0, Lj>Q'^, 

enable us to rewrite Eq. ( l25l) as 

1 f^? 

«m(i) = - / dui r]t,uJiCoth {uji,h[3/2) cos [uJe,{t)] , (28a) 
Jo 

1 r^^f 

"n(i) = / do;^ r^^tu^sin [a;£(t)] 

= ^(nFco.[n?w]-5^!M). (28b) 

As we are interested in times much longer than the typical value 1/^^^, it follows that 
sin (t)] /nt ^ 5 (t) and, consequently 



* sin [ri?(r 



d r _ y^l^r)] ^ 0, (29a) 







TTT 



,tl.,m^yM ""[f'^V" - (29b) 

TT (r — rj 2 
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With these approximations, the propagator becomes 



J {{^e} , {Vi} ^^A^'i} Ay'e} ^'^) 



exp 



2 f2 2 I /2\ 

Xfi He I Ve ) 



D{x,}D{y,} 



X exp < — 
n 



i^sA^^^ ~ ^sAye]) + Ss^/S2 [{xi}] - Ss,/S2 [{ye}] 



X exp <^ - ^ 



h 



^0 



- dr dt' [xe{T)-ye{T)][xe{t')-ye{t')] 



X / due uj£Coth.{uiihPi/2) cos[uJi{T — t')] 



~i dr [xe{T)ye{T) - ye{T)xe{T)] 



(30) 



where we have defined the relaxation constant 7^ = rif/2^u for the two system- reservoir 
couphngs and the renormahzed actions 



SsA^e\ = dr Ci{xi,xe,t) 



dr 



- Ve{xe) 



(31a) 



^s^/s-ili^i}] = / dr Li2{{xi} ,{xe} ,t') 
Jo 

= - dr [A11X1X2 + fii2 {xi - \21X2) {x2 - \12X1] 
Jo 

+Ai2/U22a;i (x2 - ""^^i ] + -^2iAiiia;2 (ii - -^^2 



(31b) 



with similar expressions for the variables ye. The renormalized potentials Ve{xe) = 
me [ujj — (Acijf)^] x'j/2, follow from the system-reservoir couplings which induce the shifts 
(AuJe)'^ = 2rieQf /Trnie. Thus, we can define a renormalized frequency, given by ujj = 
coj - (Acoefm. 

From the result that the functional integral for an infinitesimal time evolution can be 



approximated by 



33|, 



Dgexpl 15)^1 exp Us 



(32) 



we consider the evolution of the reduced density operator p on the infinitesimal time interval 
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between t and t + e (e — 0), proceeding from Eq. (fT9|) . to obtain 

pi{xe} AVi} + 

d {x',} d {y',} J ({x,} , {yj ,t + 6- {x',} , {y',} , t) p{{x',} , {y',} , t) . (33) 

We also assume the high-temperature hmit ksTi ^ hwi (for both reservoir frequencies 
uj£ <^ ^f), which allows analytical solutions of the integrals (in the variable ui) when 
defining the propagator J and using, for any function F, the approximations Xe ~ Pu/£, 
ye ~ P2e/£, and 



/t+e 
dTF[xe{T)] ^eF 



Xi + x'^ 



(34) 



where we have defined the variables Pu = — and i32e = Ve — y't With the above 
approximations, we obtain from Eq. (l30l) the propagator 



J {{^t] , {Vi} ,t + e;{xi- Pu} , {ye - p2e} , t) 



M exp { - 



^ J2 (-1)' PiiPi2 + e (An + A12A21 



/^12j 



nl P2e 



e'+e" 



n I P2e\ n ( Pie 

P21" \ yi - — \ - Pu" ^ 



+-Aa//i^'j 



X JJ^exp 



h 



[Pu - Pie) 7. 



It {xe - ye) {Pie + p2e) 



le 



T/ / P^^\ r/ f P"^^ 

Ve{xe - — j - Veiye - — 



P 



_ . le 
2 



- y. 



n 



Xi 



p 



u 



Pn 

2 



n 2 



(35) 



where A/" is a normalization factor. Since the fast-oscillating terms in the integrals in Eq. 
[32]) contribute only for P^k ~ a/^V^? '^^ expand both sides of Eq. ( !32l) up to terms 0{e) 
2]. Proceeding to a further change of variables: /?^^/ = Pu> + (— 1)^ 7£'(a;£' — ye')^-, keeping 
again terms up to 0{e), we obtain for the zeroth order term, the normalization factor 
M = (27rfe:)^ / (/i^2 ~ /^i/^2), and for the first order term the desired equation of motion ^ 
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dp 
di 



92 



+ le [xe - ye) 



2me \dx'j dyf J ' ^ " \dxi dyt 



d d 



+ 5: 



92 52 



xe - ye) 



A I \ .A^'^m^7£ X .All , N 



where we have defined the effective damping rates 

7^ 



le 



(36) 



(37) 



1 — A22'^1^2 ' 

which increase with increasing couphng strength A22- The operator equation associated with 
the above coordinate representation turns out to be 



dp 
di 



-|[H„p] + ^7, 



if f ^1 2mikBTi 

- [xe, iPe, p}\ ^ [xi, [xe, p\\ 



(38) 



5Z (-^22 [xe, {Pe',p}] + [xe, {xe'.p}]) 



h 



i'i^e) 



where the Hamiltonian Hj, is given by 



Hj = i ^ ( — + meujxU + A11X1X2 + X22P1P2 + A2iPia;2 + Ai2Xip2- 
2^\me J 



(39) 



Note that, when turning off the couphng between the two oscillators, we obtain from 
Eq. (1381) two independent dissipative oscillators described by two independent CL models. 
Therefore, in the case of distinct reservoirs, there is no effective coupling induced between 
the oscillators. 

B. A Common Reservoir 



For the case of a common reservoir, we obtain the influence functional 
J^[{xe},{ye}] = J dXdX'dY'pn (X',Y',0) J DX(r)DY(r) 



X exp 



h 



Sr [X] - Sr [Y] + {^s,/R [xe, X] - Ss,/R [ye, Y] 



(40) 
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which differs from that of the CL model for an additional system-reservoir coupling. As- 
suming identical system- reservoir couplings Cik = Cok = Ck, the solution of the influence 
functional (1^ . following directly from that in Refs. |33|, [SJ], is given by 



J" [{Xi} , {Vi}] = exp (^^Ylj^ dr^ dt' [x,{t) - yt{T)]aK{T - t')[xi{t') - y,{t')]^ 
X exp f--^^^ dr^ dt' [x,{r)-y,{T)]ai{r-t')[x,{f)+yi{t')]y 



(41) 



(^2 

a^{t) = V ^ coth {ujh(3/2) cos [uj{t)] , (42a) 
^rrikUJk 

k 

ai(t) = -5^-^sin[^(t)]. (42b) 
Defining, as in the case of distinct reservoirs, a spectral function 

Xico) = Mco) = J2 l^K^ - ^k) (43) 

constrained to the range defined by a frequency cutoff fi*^, and also a damping constant rj, 
we obtain for an Ohmic reservoir 

Xiu^) = { ' ' (44) 

0, u> Q^. 



Using exactly the same approximations performed for the case of distinct reservoirs, we 
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obtain the influence functional 



^ 77f Jo 



— exp 



T)xe>{T) - ye{r)ye'{T)) 

+l^-f f dr lY{xe{r)ye>{r) -ye{r)xe{r)) + ^ {xe{r)xeir) - yeir)ye'ir)) 
f dr !\t' [x,{T)-y,{T)][x,{t')-y,{t')] 

Of, JO JO 



X exp 



c 



X / duj ijJcoi\i{ujhj3/2)cos[ijj{T -t') 
'o 



where 



F ({x,} , {x',} , {y^} , {y',}) = H^"^ ' ' W ' ^S] 



(45) 



(46) 



Wc stress that wc have defined, after the assumption {Cik} = Ck, the relaxation constant 
7 = rj/2ij, for both system-reservoir couplings, which implies immediately that {/J^ee} = iJ- 
and, consequently, {m^} = m. For the master equation, we obtain 



.h\22 



E 



2m'ykBT 



{- 


--)] 


\dxe 


dyij . 



d d 
dxe oyi 



2 \ dx^dxei dyedy^ 
+^Aii [xtxv - ytyv) + 1 — ^ — [xt + y{) ^ [xe' - ye") 



h 



(47) 



in which the renormalized damping constant 7 = 7/(1 — m^A|2) and the effective coupling 
parameter All = An — 2QPri/T: are considered. As usual, this shift for the effective couphng 
parameter is small and can be included in An. The operator equation associated with the 
above coordinate representation turns out to be 



dt 



h 



-E 



LV 



■.7 r r -11 , 2m-fkBT ■ 
^^ R, {Pe':p}\ + ^ R, [xe', p\\ 



(48) 
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where the Hamiltonian Hjj, given by 

= ^ ( — + mwlxf] + A11X1X2 + A22P1P2 + A12X1P2 + A21P1X2, (49) 
contains the renormahzed frequencies 

wl = u)l + 27 + A2i5^2) • (50) 

A reservoir-induced coupling between the oscillators 

It is worth stressing that, differently from the case of distinct reservoirs presented in Eq. 
( 137|) . the common reservoir induces an effective coupling between the two oscillators, even 
when their original interactions {Xw} are turned off. In fact, with X^' = 0, the Eq. 
simplifies to 



fj fll 



[x£, {P£',p}\ H ^2 F^' F'?'' PJJ 



(51) 



giving two independent CL models, for i = i', apart from the reservoir-induced coupling 
between the oscillators, for i 7^ i'. This effective coupling, for the case analyzed here of 
an ohmic reservoir in the high-temperature regime, consists therefore of a dissipative and 
a diffusive term, given by (j/ih) ([xi, {p2,p}] + [x2, {pi,p}]) and {2m'~^kBT /h^) [xi, [xi',p\], 
respectively. Interestingly, such dissipative and diffusive terms couple together the variables 
of both oscillators while, evidently, the equivalent terms in the CL model apply to the 
variables of a single particle. 



V. SOLUTIONS OF THE MASTER EQUATIONS dSHD AND (HT]) 

In this section we present the solutions of the master equations governing the dynamics 
of the coupled dissipative harmonic oscillators in both cases of separate reservoirs and a 
common one. These solutions enable us to analyze the coherence and decoherence dynamics 
of quantum superpositions prepared in one of the oscillators of entangled states prepared 
in both oscillators. Moreover, the form of the solutions presented here enable a complete 



understanding of the evolution of such states, thus en 
and decoherence analysis offered by the CL model {4]. 



arging the perspective of the coherence 
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A. Distinct Reservoirs 



The introduction of the collective (Re) and relative (re) coordinates 

Xi + yi 



Re 



re = Xi- ye 



(52) 



enable us to rewrite the master equation (!36|) in the form 

dp 
di 



~ ) (-1) ^— Q,. an - 27^^"^— + i^—Rere - —r^ 



mi dr^dRi 



' dre h 



-E 



oRi' or I I oriORi' n 



(53) 



where the diffusion coefficients Di and the effective coupling parameters Ti and are given 
by 



Di = 2mi-fikBTi, 

= All + 2mi^i ^ Xi'i, 

Ai = 2X22mi^i + ^ Xii'. 



By the partial Fourier transform 



1 / ;>+oo 



PU^^M^^I,IJ = — I I I / ai?,e-*^^^^ lp({i?a,{r,},t) (55) 
we reduce the second order partial differential equation (!53|) to the first order one 



(54a) 
(54b) 

(54c) 



dt I ( 7f ^ rj^^ + 



d iTLiuj d Di 2 



9r/? h dKf h? 



- Aa//^^/ ) — — + {AifTi' - hX22Kp) — 
n J oKf ore 



whose solution can be obtained by the method of characteristics 
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(56) 



36| . Defining the curves 



Ki = Ki{s), Ti = ri{s), and t = t{s) 



(57) 
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we obtain, from the partial differential equation (156|) . a system of coupled ordinary differen- 
tial equations 



ds 



A22 (1 — Sw) H 6u' 

nit 



d_K^ = ^^^^^^m + ^ [(Pars - hXi2K2) + (Firi - h\2iKi) 81,2] 
as n h 
dp 

ds " 



ds 

The first four equations of the above system can be expressed, in matrix form, as 



(58a) 

(58b) 
(58c) 

(58d) 





1 n [t) \ 


/ 


271 


h_ 

mi 


A2 


— hX22 


\ 




d 


Ki (t) 









^2 






Ki (t) 


dt 


r2 (t) 




Ai 


— h\22 


272 


h_ 




r2 (t) 




\ K2 (t) 1 


\ 


h 


— ^21 


\m2Ujl 





/ 


V K2 (t) ) 



(59) 



where we shall denote the square matrix by M and its eigenvalues by Am- Next, we note 
that the solution of this system of coupled equations can be written in the form 



/ n it) \ 




' Vn V12 Vi3 Vi4 ^ 




( ci it) \ 


Kiit) 




1121 V22 V23 V24 




C2 it) 


r2 it) 




VSI V32 V33 V34: 




C3{t) 


\K2{t)) 




\ V41 V42 Vi3 Vii ) 




\ C4 it) j 



(60) 



where the eigenvector ^ rji^ i]2n ^sn V4.n ) associated with the eigenvalue A„. 

Making use of both solutions for (t) and that for the third ordinary differential equation 
fl58cp . given by 



p{{Ke},{r,},t) = Bexp 
we finally obtain the density matrix 



(61) 



~p{{K,} , {rj , t) = p{{K',} , {r;} , 0) e-^d^^.W},*)//^^ , 



(62) 
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with 

Z ({ira , , t) = 5^ C„ it) C„ (t) Dr^^rnmn + D.r,,„^V.n _ ^-(A„+A„),^ ^ ^gg) 

and Cm if) = Cm (0) e^™*, where (0) is determined by the initial conditions. 

The next step is to calculate the inverse Fourier transform of Eq. ( I55l) . To this end, 
let us suppose that the coupled harmonic oscillators are prepared at t = in the general 
superposition of Gaussian functions of width {a^}: 

m ({x J , 0) = ^ JJ ^ (64) 

which include both separable and entangled states, depending on the choice of the Gaussian 
centers qu and q2e for oscillators 1 and 2. Rewriting the density matrix for this wave function, 
given by 

Pi{xe} , {ye} ,0) = J2 P^^t' n e-[(^^"+'^^"^)/"^"]' e-[(^^"+«^"^)/"^"]' , (65) 
1,1' I" 

with the collective and relative coordinates defined in Eq. fl52l) . and computing its Fourier 
transform, we obtain 



p(ra,W},o) 

^^P^'P^" J^cr^exp { -2 



qw - qu" + r'l, 

1(7, 



. (66) 



Therefore, from the results in Eqs. (J62I) and (16611 . we obtain the general solution for the 
transformed master equation (!56|) 



p ({Xa , {rj , t) = 7^ ({i^a , {r,} , t) J]T«, 



Ll' 



X exp 
X exp <! % 



(67) 



which, for t = 0, gives the initial condition (l65l) . Assuming that the matrix elements 
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follow from the inverse of the matrix composed by r].nm, we have defined the functions 
n {{Ke} , {re},t) = exp {- [$ii(t)r? + $33 ^rs' + $i3(t)rir2 + ^laWrii^i + ^uit)nK2 

+ ^22{t)Kf + $44(^)^-2 + $24(t)Klir2 + ^23{t)r2Ki + ^u{t)r2K2] } 

(68) 



T«/ = ^PePe' Yl ^i" exp ~( 



a III 



(69a) 



id=Uj>i) 



+20, (1 



-(Ai+A2)t 



)] 



£ik£jk' — (1 — <^A:fc') f^jfc'^jfc] 



Q'w - qw . q2e - q2i' \ _A,t 
2 Vij + -2 Vsj ) £jk e 



(To 



qu + qw , q2e + q2i' \ -A^t 
n + ^ V4j £jk e 



(69b) 
(69c) 

(69d) 



The next and final step is to calculate the inverse Fourier transform of Eq. fl67j) . At 
this point we observe that it is highly recommendable to write the Fourier transformed 
density matrix p ({i?^} , {r^} , 0) in its initial normalized form, obtained by a coordinate 
transformation of Eq. fl65|) : 

-1 



p {{Ri} , {rj , 0) = <^ 7ro-i(T2 



exp 



X exp 



1 / qei - qe2 

2 V ff£ 

n2 



jqee' - qu" + re) + {qw + qu" + ^ReY 
2a2 



(70) 



The exponentials in Eq. ( ITOi) associated with i = i', represent the diagonal elements of 
p{{Re} , {ri} , 0), which are associated with the probability amplitudes of the system, while 
exponentials associated with i ^ i' represent the off-diagonal elements, associated with phase 
coherence. After the sudden system-reservoir couplings it is expected that the diagonal 
elements will be dragged to the origin of the coordinates by the dissipative mechanisms, 
whereas the off-diagonal elements will vanish continuously, owing to the associated noise 
injection into the system. Therefore, by keeping exactly the form of the above initial density 
matrix (170|) after its time evolution, we directly verify such expected dynamics, simplifying 
our evaluation of the decoherence effects. From this perspective, we obtain the final solution 



22 



1+n 



exp 



1 / qn - q£2 

2 I ai 



-1 



X exp 



--Eee {{Re} , {re} , t) + iQee {{Re} , {re} , t) 



(71) 



where the functions E^n{{Re} , {re} ,t), accounting for the dynamics of the diagonal and 
off-diagonal elements, are given by 

-mn{{Re},{re},t)-y^2u 

,5„,-ire - [2 {qei{t) - %2(t))]'^^' re' - [2 {qen{t) - qe'2{m'"' 



+2" 



(^ee'{t) 



(^e'e{t) 



, (72) 



whereas the functions 

e^. {{Re} , {re} , t) = J] << 2 [i?, + CH^)] 



-2)' 



^12(i)ri+$23(0'^2 + ^^U0 



+ (-2)' 



$i4(t)ri + $34(t)r2 + 9il{t) 



T.l{t) 
+ O^Z-'\t)re). 



(73) 



account only for an oscillatory dynamics. Moreover, we also have, as part of the normaliza- 
tion factor 



E(t) = ^4$44(t)*22(t)-$i4W, (74) 

the functions associated with the widths of the diagonal and off-diagonal Gaussian peaks: 



23 



Sii(t) 

S22(t) 



/2S2(t) 

$44 (t) 

/2S2(t) 

$22 (t) 



Sl2(t) = S2l(t) 



/2S2(t) 

$24 (t) 



(Tn(t) = W2$n(t) -2 



$22(t)$?4(^) + $44(t)$?2W - $12(t)$14(t)$24(t) 



S2(t) 



f^22(t) = W2$33(t) -2 



$22(t)$i4(i) + $44(t)$i3W - $23W$34(t)$24(t) 



0-12 (t) = 0-2l(t) 



2$13(t) - 2 



2$22(t)$14(t) - $12(t)$24(t) 

S2(t) 



$34 (t) 



^ 2$44(t)$12(t)-$14(t)$24(t) ^ 



(75a) 
(75b) 
(75c) 
(75d) 
(75e) 



(75f) 



Finally, regarding the exponential decay multiplying the off-diagonal elements of the 
reduced density matrix (I7T1) . we obtain the function 



m 

where 



E 



qn — Qi2 



E 



'it)m +5:2^--^ (-1)^-^' ;f " 



(76) 



(1) [2$44(t)$12(t) - $14(t)$24(t)] ef^{t) + [2$22(t)$14(t) " $12 (t) $24 (t)] ^iJ (t) 



S2(t) 



^2(t)=^^i^(t) 



(77a) 

(3) [2$44(t)$23(t) - $34(t)$24(t)] ^2^) + [2$22 (t) $34 (t) " $23 (t) $24 (t)] (t) 



(77b) 



Comparing the initial and the evolved density matrices, given by Eqs. ( 1701) and (|7T]) . we 
verify, as expected, the displacement of both the diagonal and off-diagonal peaks towards 
the origin of the coordinates, described by Eq. ( |72ll . Apart from this displacement towards 
the origin, the system-reservoir coupling also induces the oscillatory term defined by Eq. 

m. 
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A general decay function of the off-diagonal peaks of the density matrix 

From the evolved density matrix in Eq. flTTl) . we also deduce a general expression for 
the decay 'D(t) of the off-diagonal peaks as time goes on, represented by the exponential 
function 

V{t) = exp |-2 [1 - r(t)] ^ [{q,^ - fc) /2o-,]'| . (78) 

This decay function offers complete information relative to the decoherence of any initial 
state of the coupled dissipative oscillators. We anticipate that the same form of the decay 
function in Eq. (1781) will apply to the case of a common reservoir, except for the time- 
dependent function T{t), which differs in the two CcLSGS cLS discussed below. 

B. A Common Reservoir 

The solution of the master equation for the case of a common reservoir is entirely analo- 
gous to that of separate reservoirs. From the early assumption of identical system-reservoir 
couplings Cik = = Cki rendering the same masses m, damping rates 7, and diffusion 
coefficients D for both oscillators, the master equation, rewritten in terms of the collective 
and relative coordinates in Eq. ( l52l) . is given by 



dp 
di 



h 92 



D 



m orfoRf ore n 



drfdRf h 



z—reRe 



(79) 



where the diffusion coefficients Di and the effective coupling parameters and are given 
by 



D = 2m'ykBT, 

= All + 2m7 ^ Xe>i, 

Ai = (X2iSi2 + XuSn) + 2m7A22- 



(80a) 
(80b) 

(80c) 



Following the steps outlined above for the case of distinct reservoirs, leading to the system 
of coupled ordinary differential equations ( IHTi) . we now obtain the analogous system 
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ds 



A22 (1 - ) H 5. 



m 



dKi mujjri 1 



ds 



+ ^ [(r2r2 - hXi2K2) 5n + (Firi - ^A2i-ft'i) fe] , 



ds 

As before, the first four equations is given in the matrix form as 



(81a) 

(81b) 
(81c) 
(81d) 





1 n it) \ 


1 


27 


_A 

m 


A2 


-h\22 


\ 


/ n (t) \ 


d 


Ki (t) 











— A12 




^1 it) 


dt 


r2 (t) 




Ai 


-hX22 


27 


_A 

m 




r2 it) 




\K2it)J 






— A21 







/ 


\ K2 (t) 1 



(82) 



From here on, all the derivations performed earlier for the case of distinct reservoirs can 
be followed in exactly the same way, leading to the same structure as the solution given 
in Eq. (|7T1) . The difference is that the elements of the above square matrix lead to values 
of Tjmn that differ from those in the case of distinct reservoirs. Moreover, the function 
Z {{Ki} , {r^} , t) becomes 

Z {{Ki} , {n} ,t)=DY,C^ {t) Cn (t) (^1- + (^m + VSn) _ ^^(A„+A„)t) ^ (gg) 

m,n 

instead of the expression given in Eq. fl63|) . Therefore, in spite of the solutions to the master 
equation having the same form, the evolution of the reduced density matrices of the two 
dissipative oscillators must be quite different for the cases of one common reservoir and two 
separate ones. We finally point out that the expression for the decay or decoherence 'D{t) 
of the off-diagonal peaks of the initial density operator also has the same structure as that 
in Eq. fl78|) . the only difference being the time- dependent function T{t), which differs in the 
two cases of a common and two separate reservoirs, again due to the values of rjmn- 



VI. DECOHERENCE IN THE DOUBLE CALDEIRA-LEGGETT MODEL 

In this section we analyze the decoherence of three particular entangled states prepared 
in both oscillators of the network at t = 0. These entanglements, derived from Eq. (IM]) 
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under the assumption that they all have the same initial mean energy {E) = {Hs^ + -^^^2) 
and distance d between the two peaks in the positional a;iX2-space, are given by 



{{xe} , 0) = Ar(i) jexp - (x, + ^) ' /a? 



+ exp 

^(2) , 0) = Ar(2) <; exp 

+ exp 

^(3) {{xe} , 0) = AT^^) <; exp 



exp 

2 



Xi + Qi ) 



exp 
2 



+ exp 



exp 



3^2 -g?^) /(rf 

(2)V/ 



(84a) 



(84b) 



(84c) 



The choice of equal mean energies (E) and distances d follows from the fact that the deco- 
herence process depends on the energy, the distances between the components of the super- 
position, and the damping rate defining the system-reservoir coupling. It is well-known that 
the decoherence time varies inversely with the energy, the distance d, and the damping rate. 
The above states differ from each other only by the position of their peaks in the XiX2-space. 
While both peaks of {{^e} ? 0) (\E'(2) ({a;^} ^ 0)) are positioned in the fourth (first) quad- 
rant of the XiX2-space, those of the state "^^^^ {{xe} , 0) are positioned in the first and third 
quadrant. Moreover, whereas "^^"^^ {{^i} , 0) is obtained from {{xi} , 0) by a rotation in 
xiX2-space, {{xi} , 0) requires, apart from the rotation, a displacement operation over 

We analyze the decoherence time of the above entanglements for the two cases of a com- 
mon and separate reservoirs. We observe, for comparison, that the analysis of decoherence 



in Refs. 



10 



13 



14l |. centered on absolute zero reservoirs, shows that the decoherence rate 



for the case of a common reservoir is attenuated, compared to the case of distinct ones. 
However, the present analysis is based on the opposite scenario of the high-ternperature 



limit, so that we do not expect to obtain similar results to those in Refs. 



10|,ll3|,ll4| 



As the decoherence rate Vit) is given by Eq. (1781) . in Figs. 1, 2, and 3, we plot Vit) 
against the scaled time 7it for the states \E'*^^\ and respectively. In parts (a) and 
(b) of Figs. 1 and 2 we plot the decoherence rate T>{t) for the cases of distinct reservoirs 
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and a common one, respectively. Adopting unit constants h,kB = 1, masses (mi = 1712 = 1), 
frequency uji = 1 and widths ui = (12 = 1, we have set the magnitudes Ti = T2 = lO^hwi/kB 
and q[^^ = q!^^ = 2q!^^ = 2q^^ = lOcri. We have also assumed a regime of parameters where 
7i = 72 ^ {^u'} < uji = UJ2/2, with 71/^011 = 10~^. With these values, the mean energy 
(E) and distance d associated with states fl84p becomes (E) ^ 158huJi and d = ^^plox. To 
reach the same (i?) and d for all three states we have assumed that a\ = 12as, apart from 
the relation gf ^ = ^/3/2ql^^ = ^/l5/2ai. 

For comparison, the thick solid line in all three figures represents the decoherence time 
of the Schrodinger-cat-like state 

* (x, 0) = {exp [- (x + qf /a^] + exp [- {x - qf /a^] } , (85) 

prepared in one of the oscillators, decoupled from the other. This Schrodinger-cat-like state 
also leads to the same values established above for the mean energy (E) and distance d 
between the two peaks in the x-space. To achieve this, we set the relations q = 5cri/v^and 
cr = 6 X lO^^fTi. Therefore, the thick solid line describes the decoherence process of a 
superposition state in the CL problem. 

In all three figures, the solid (dashed) and dashed-dotted (dotted) lines describe the 
decoherence processes when considering distinct reservoirs (a common one) and the coupling 
between the oscillators given by Aiigig2 + A22P1P2 and Xi2qiP2 + A2i5'2Pi, respectively, with 
Xw = 0.1. Both couplings Xuqiq2 + X22P1P2 and Xi2qiP2 + X2iq2Pi, when described in terms 
of the usual annihilation (creation) operators 01,02 {a\,al), correspond to the rotating terms 
a\a2 + (21(22 and the counter-rotating terms i (^aia2 — a\al^, respectively. We observe that 
the decay rates of the curves in Figs. 1 and 2 (a and b) are around that associated with the 
Schrodinger-cat-like state for both cases of a common and distinct reservoirs. Therefore, the 
case of a common reservoir does not exhibit advantages over that of distinct reservoirs, as 
de.on...ed p.™u.y fo. absolute .e.o .e.™. BflQ. 

Considering now Fig. 3, we observe that all curves decay faster than those in Figs. 1 and 
2. Moreover, the decay rates of the curves associated with the coupling Aiigig2 + X22P1P2 
are even faster than those associated with Xi2qiP2 + X2iq2Pi- This behavior follows from Eq. 
(1371) . which shows that the larger the coupling strength A22, the larger the effective damping 
rate 7^. The same explanation applies to Fig. 4, where the same curves as in Fig. 1 are 
plotted for larger strengths Xw = 0.5; we observe that the decay rates of the curves derived 
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from the coupling Aiigig2 + are significantly faster than those for XuCliP2 + ^2192^1, 

as in Fig. 3. In Figs. 1 and 2, the effect coming from the relation between A22 and 7^ is 
blurred, making the decay rates of the curves associated with the coupling \11q1q2 + A22P1P2 
similar to those for Xi2(liP2 + X2i(l2Pi- 

Differently from the works in Refs. 10, [isl, Q], where the decay rates were chiefly de- 
termined by the assumption of one common or two separate reservoirs, here the difference 
between the decay rates comes from the different coupling mechanisms between the oscilla- 
tors, Aiigig2 + A22P1P2 or Xi2qiP2 + A2ig2Pi- The reason for this is that Refs. [iol,Q,Q 
apply to absolute zero reservoirs, where the coupling between the oscillators induced by a 
common reservoir tends to decrease the decoherence rates. The same effect of coherence 
control is also achieved in Refs. 



10|,[l3|,Q dilators are 

strongly coupled to each other. Especially in Ref. it is demonstrated that the R-DFSs 
emerge from situations where the whole network interacts with a common reservoir or when 
each resonator, strongly coupled to each other, interacts with its own reservoir. The present 
work, however, applies to the high-temperature regime where, as demonstrated in Figs. 1, 
2 and 3, the interaction induced by a common reservoir is completely blurred by the high- 
temperature effects. From this fact we may expect the high-temperature regime to prevent 
the emergence of R-DFSs, at least in a network with a small number of oscillators, as in the 
case at hand. 



VII. CONCLUDING REMARKS 

In this study we analyze the double CL model, i.e., the path integral approach to two 
interacting dissipative harmonic oscillators. We derive and solve the master equations asso- 
ciated with two different situations: when each oscillator is coupled to its own reservoir, and 
when both oscillators are coupled to a common reservoir. In both cases, the derived master 
equations consist of two CL models — describing the dissipative mechanism of each oscilla- 
tor independently — as well as the dynamics arising from the interaction between the two 
oscillators. However, in the case of a common reservoir, we identify a reservoir-induced cou- 
pling between the oscillators, even when the original interaction between them is schwitched 
off. Such a reservoir- induced coupling, recently pointed out in Ref. [28|, encompasses both 
dissipative and diffusive terms which couple together the variables of both oscillators. These 
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terms thus account for the energy loss of the oscillators through each other, apart from a 
joint diffusive process. 

The occurrence of such a reservoir-induced coupling between the oscillators was also 



3, Q, Q, ^ 

oscillators are treated through perturbative master equations. In Refs. [id, llSl, l2Q|] , the 



pointed out in Refs. 

m lisl, y, l2o|, where networks of dissipative quantum harmonic 



e quantum n; 

i. Q,Q,QJ 



occurrence of indirect-decay channels — by which the network oscillators lose excitation 
through all the other oscillators — is demonstrated in two different situations: when all the 
(non-interacting or interacting) oscillators are coupled to a common reservoir and even when 
strongly interacting oscillators are coupled to their own reservoirs. 

Regarding the master equation solutions reached in this paper, we stress that their form 
enables a full comprehension of the evolution of initial states of the network, enlarging 
the perspective of the coherence and decoherence analysis offered by the CL model {4]. 
Through these solutions we compute a general expression for the decay rate of the off- 
diagonal peaks of the density matrix of initial superposition states, which applies to both 
cases of a common and distinct reservoirs. Such an expression offers complete information 
about the decoherence of the initial state of the coupled dissipative oscillators. In this regard, 
we have analyzed, as an application, the decoherence process of particular entanglements 
in the positional space of both oscillators. The results demonstrated that the coupling 
induced by the common reservoir does not lead to the collective damping effect mentioned 
above. The high-temperature regime of validity for our calculations completely blurred such 
reservoir-induced coupling which at absolute zero works, in general, to delay the decoherence 
process or even to produce the R-DFSs. However, we find that different interactions between 
the dissipative oscillators, described by rotating or counter-rotating terms, result in different 
decay rates of the interference terms of the density matrix. The decay rates associated with 
the counter-rotating terms of the interaction between the oscillators are significantly faster 
than those coming from the rotating terms. The reason for this is that the effective damping 
constants increase with increasing coupling strength related to the counter-rotating terms. 

We note that a recent paper addressed the question of the derivation of a master equation 
for two coupled harmonic oscillators through the influence-functional method of Feynman 
and Vernon js?]. However, the authors assume a coupling between the oscillators different 
from that given by Eq. (jl]), apart from considering only the case of a common reservoir. We 
finally stress that the present double CL model can be used for the analysis of dissipative 
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bipolarons, besides other problems in several areas of physics where the CL model has 
been successfully employed over the last few decades. A more detailed analysis of the 
development in Ref. [28| is one immediate application of the present study. The effect of 
temperature on decoherence and the emergence of R-DFSs in the low-temperature regime is 
also a point worth looking into. Furthermore, the tunneling process of coupled dissipative 
systems is a phenomenon which may be accounted for by the present work. Together with 



the achievements in Refs. 



10 



13 



14 



20| , we believe that the present work furnishes a great 



deal of material for discussion of the physics of coupled dissipative systems. 



VIII. APPENDIX A - DIAGONALIZATION OF THE HAMILTONIAN i^Si+Sa 



The Hamiltonian modeling the two harmonic oscillators and their mutual interaction, 
Hsi+S2: &s given by Eq. (IT^ with Vi{qi) = miUjjqj/2, can be rewritten in terms of the 



quantum mechanical creation al and annihilation ae operators, as 



Hsi+S2 — ^ 
with the coupling strength 



uJe (^alac + ^ + {giaia2 + g2aial + H.cj 



(86) 



All - 2A2imiCc'i , / , , . Ai2 \ 

ge = — - + (-1) y/miUJim2UJ2 X22 + i • 

y/miUJim2UJ2 V rriiUJi J 

The diagonalization of the form (!86|) . easily performed than that in Eq. 
Hsi+S2 = f>'J2e^i (^^l^i + i)' with the normal-mode frequencies. 



(87) 

(fT2D, leads to 



^2 _ ^1 + ^2 I _ |2 , I _ |2 



l^ll 



-(-1 




and the normal-mode operators 



,2\ 2 



l^fil^ {uJi - uj2f + 192^ (t^i + ^^2)^ (88) 



(89) 
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where 



X 



(90a) 



(90b) 



l^ll' - 1^2!' + {^l" - ^1) {^i" + (-l)'+''c^2) 

The normal-mode coordinates and P^, following from the operators and ^4^, are given 

by 



'2n 



with the coefficients 



2H 



and phase factors (f)e = [Oi — (— 1)^6'2] /2, where the quantities 

9(, — cos~^ 



(91a) 
(91b) 

(92a) 
(92b) 



U + {-l)%2 



[61 + (-1)^62]' [62 +61 

are defined by the dimensionless strengths 



2 r ' 



(93) 



— 



{mim2ijJiUJ2f^'^ 5a + 5, 



i2 



(94a) 



2 Y m£iuj£' 



(94b) 



Through the coordinates Qe and P^, we finally obtain the diagonalized Hamiltonian 



(95) 
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Figure captions 

Fig. 1 Plot of the decay or decoherence function 'D{t) against the scaled time 7it for 
the states in Eq. fl84al) . considering the cases of (a) distinct reservoirs and (b) a 
common one. The thick solid line describes the decoherence process of the Schrodinger- 
cat-like state \l/ (x, 0) in Eq. fl85|l . prepared in a single dissipative oscillator. The solid 
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(dashed) and dashed-dotted (dotted) lines describe the decoherence processes of when 
assuming distinct reservoirs (a common one) and the couphng between the oscillators given 
by Aiigig2 + ^22PiP2 and \12q1P2 + ^2152^1, respectively, with A^/ = 0.1. 

Fig. 2 Plot of the decay or decoherence function 'D{t) against the scaled time 7it for 
the states in Eq. (184bp . considering the cases of (a) distinct reservoirs and (b) a 
common one. The thick solid line describes the decoherence process of the Schrodinger- 
cat-like state \l/ (x, 0) in Eq. (153]) . prepared in a single dissipative oscillator. The solid 
(dashed) and dashed-dotted (dotted) lines describe the decoherence processes of when 
assuming distinct reservoirs (a common one) and the coupling between the oscillators given 
by Aiigig2 + A22P1P2 and Ai2giP2 + A2ig2Pi, respectively, with A^/ = 0.1. 

Fig. 3 Plot of the decay or decoherence function 'D{t) against the scaled time 7it for the 
states ^^^^ in Eq. fl84cp . considering the cases of distinct reservoirs and a common one. The 
thick solid line describes the decoherence process of the Schrodinger-cat-like state \1/ (x, 0) in 
Eq. (|85l) . prepared in a single dissipative oscillator. The solid (dashed) and dashed-dotted 
(dotted) lines describe the decoherence processes of when assuming distinct reservoirs 
(a common one) and the coupling between the oscillators given by Aiigig2 + A22P1P2 and 
Ai2?iP2 + A2iq'2Pi, respectively, with A«/ = 0.1. 

Fig. 4 Plot of the decay or decoherence function T>{t) against the scaled time 71/!: for the 
states in Eq. fl84al) . considering the cases of (a) distinct and (b) a common reservoir. The 
thick solid line describes the decoherence process of the Schrodinger-cat-like state \1/ (x, 0) in 
Eq. fl85l) . prepared in a single dissipative oscillator. The solid (dashed) and dashed-dotted 
(dotted) lines describe the decoherence processes of when assuming distinct reservoirs 
(a common one) and the coupling between the oscillators given by Aiigig2 + A22P1P2 and 
Ai2Q'iP2 + A2iq'2Pi, respectively, with A«/ = 0.5. 
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